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TECHNICAL MEMORANDUM X-53767 


FLOW PARAMETERS OF THE CORONA 
SUMMARY 


The coronal plasma is assumed to be a fully ionized mixture of two species. 
Nondimens ional parameters are obtained from the fundamental equations of plasma 
flow . Some thermodynamic properties of the coronal plasma are presented in 
tables which are based on a coronal model computed by F. L. Scarf and L. M. 
Noble. 


INTRODUCTION 


A steady-state, spherically symmetric fluid-dynamic model under quiet- 
day solar conditions is adopted for the corona. A steady-state quiet solar wind 
blows continuously with the streaming becoming supersonic at great distances. 

The local heating energy is transferred primarily by thermal conduction. The 
corona is assumed to be a fully ionized mixture of two species. Nondimensional 
parameters are obtained from the following fundamental equations, which describe 
the plasma flow: 


P = — PT 
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e e 
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In the equations for the coronal plasma, which were^ given, the variables 
to be investigated are u,p,p , p , T , H , p or n ,J or 6 . Of these 

6 G G 

u , T , and n g are known from the adopted model. The others are unknown and 

are derived in this report. In addition some characteristic velocities, frequen- 
cies, lengths, and some nondimens ional parameters are computed. 

The coronal models of the following authors were investigated: Waldmeier, 
Baumbach, Newkirk, Whang, Liu and Chang, Parker, Chapman, Pottasch, Scarf 
and Noble. Of all models, the one by Scarf and Noble is most complete and was 
therefore adopted. The computations of the coronal parameters could not have 
been carried out with the completeness in which they are presented, using any of 
the other models. 

In the equation for F , the first term represents the solar gravitational 
field, the second electromagnetic, the third centrifugal, the fourth Coriolis, 
and the fifth external heating forces . 

In the region between 2 and 15 , where reliable density measurements 

are available, the solar corona can be described by the Navier-Stokes equations 
with no external heat source. The coronal flow appears to be stationary, radial, 
spherically symmetric, and unaffected by electromagnetic forces. The best fit 
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to n (r) is obtained for a 90 percent hydrogen, 10 percent helium corona. 

C 

Numerical values for the solution are contained in Table I. 

The streaming speed increases rapidly beyond the sonic transition near 
7 Rq . The values are not reliable outside the region 2 < R^ < 15 . 

The properties of the coronal plasma will be presented in tables . 


TABLE I. MODEL FOR A 90% H 2 -10% He CORONA 


r 

r ® 

T 

°K 

10 6 

n 

e 

m -3 

10 6 

u 

msec -1 

10 3 

1.58 

1.500 

6. 21 • 10* 

9.5 

1.99 

1.470 

2. 10 • 10* 

18.0 

2.00 

1.470 

2, 03 • 10 s 

18.4 

4.50 

1.030 

1.02 • 10 s 

72.4 

4.67 

1.020 

9. 20 • 10 4 

75.0 

6.80 

0. 889 

3. 09 • 10 4 

105.0 

7.20 

0.873 

2.62 • lO 4 

110.0 

7.25 

0. 870 

2.56 • IQ 4 

110.0 

9.56 

0. 810 

1 . 22 • 10 4 

132.0 

10.00 

0.792 

1.02 • 10 4 

146.0 

12.60 

0. 755 

4.95 • 10 3 

190.0 

13.70 

0.735 

5.00* 10 3 

161.0 

15.00 

0.735 

2 . 70 • 10 3 

246.0 

214.30 

0. 270 

6.5 

354.0 
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TYPICAL PHYSICAL QUANTITIES 


The following quantities, which characterize coronal plasma dynamics, 
can be introduced into the preceding plasma equations . 

1. Length which characterizes the dimension of the flow field. 

All distances may be expressed in terms of i ^ . For the corona, the following 
distance may be used which is larger than the scale height: 

i = 10 9 m 
o 

2. Characteristic velocity u q . The approximate velocity of expansion 
at 10 6 °K may be used as characteristic velocity: 

u =7.5* 10 4 msec -1 
o 

3. Time t which characterizes the interval considered: 

o 

SL 

t = — = 1.33 • 10 4 sec 
o u 

o 

4. Number density of the coronal plasma, n . The values are contained 
in Table II. 

5. Apparent molecular weight of the plasma, \x . With sufficient 
accuracy, it can be considered as the same as that of the Sun: 

H = 0.62 

m = 1.02908 - 10“ 27 kg 

6. Temperature of the plasma. As shown in Table I, the temperature 
varies greatly in the corona. The typical temperature may be taken with 
sufficient accuracy: 

T = 10 6 °K 
o 
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TABLE n. THERMAL PARAMETERS OF THE CORONAL PLASMA 
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7. The typical pressure p^ may differ for various flow problems . 

In table II, the densities have been obtained from the equation: 

p = mn 

and the total number density for the assumed coronal composition from the 
equation: 

n = 0.525 n 
e 

The pressures are then obtained from the ideal gas equation. The pressure 
near T q maybe considered as characteristic pressure: 

p = 2.46 • 10~ 6 kg m" 1 sec -2 
o 

8. Boltzmann’s constant k which characterizes the gas constant of the 

plasma: 

k - 1.38046 * 10" 23 joule deg -1 

kp= — = 1.35456 * 10 4 joule kg" 1 deg -1 

9. Coefficient of viscosity ju^ , which characterizes the viscous stress 
of plasma. The values given in Table II have been derived from the equation: 

5 / 

fj. =1.2* 10 -17 T /2 kg m -1 sec 1 
v to 

10. Nonelectric forces, such as the gravitational force. If the coronal 
flow is assumed to be stationary, radial, spherically symmetric, unaffected by 
electromagnetic forces, and without external heat sources, the acceleration of 
solar gravity only appears in the equation for F : 

g = 2.7398 * 10 2 m sec -2 

However, when the coronal gas possesses a sufficiently large density, its 
self -attract ion becomes important. The correct momentum equation for a sta- 
tionary, spherical, self-gravitating gas in an external field is*. 
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d 

dR 


1 2 
— mir - 
2 

GM 0 m" 

i 

d 

r 

4t r Gm p 

r 

+ 

p 

dr P 

v z J 

r 


11. Specific heat at constant volume Cy which characterizes the 
internal energy of the plasma. For a monatomic gas: 


c = — — 

V 2 m 

After substitution of the previously given values, this equation yields: 
Cy = 2 . 03 1 84 * 10 4 joule kg -1 deg -1 

12. Coefficient of diffusion D, . Due to the temperature and pressure 

d 

differences, the molecules in the coronal plasma diffuse from regions of high 
concentration to regions of low concentration, according to the equation: 

iT = - D , Vq . 
d 


In the kinetic theory, the value of is given by the expression 

ix 

D, = 1.204 — . 

d p 

The values obtained after the substitutions are contained in Table II. 

13. Coefficient of heat conductivity X, which characterizes the heat 
conduction of the coronal plasma. It is obtained from the relation: 

5 / 

X = 6 • 10" 12 T 2 joule m -1 sec" 1 deg -1 
The results are contained in Table II. 

14. Radiation constant a which characterizes both the radiation 

r 

pressure and the radiation energy. Its value is 
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= 7. 5641 • 10~ 16 joule m” 3 deg -4 . 


15. Diffusion coefficient of radiation D 

i 

radiation flux. It is derived from the relation 


. a c T* 

D = i JL^JL 

r 3 kp 


which characterizes the 


In this expression, k reduces to the electron-scattering mass absorption 
coefficient 

k = 0.019 (1 +X) = 0.0361 m 2 kg -1 . 

6 

Due to the extremely low densities, an expression of the type: 

T“ 3 . 5 

r 

K = K p 

° t 

becomes negligible. The radiation temperatures T , as evaluated from the 

effective temperature of the Sun T = 5780°K by the square-distance law, are 

© 

given in Table n . After substitution of the previously given values, the results 
contained in Table n are obtained. 


16. Magnetic induction B . At the photospheric level, the Sun has 
general induction field of 10“ 4 weber m -2 ( 1 gauss) which leads to a radial field 
of 2. 2 • 10” 9 weber m -2 ( 2. 2 • 10“ 5 gauss) at a distance of 1 AU (214 R^) . At 

the orbit of the Earth, the inclination of the lines of force to the radial direction 
becomes 56°, approximately; the azimuthal component is 3. 8 • 10“ 9 weber m“ 2 
(3. 8 • 10 -5 gauss) . Inside the orbit of the Earth, the quiet -day interplanetary 
magnetic field is predominately radial, with a tendency to spiral toward the west 
at the orbit of the Earth. Beyond the orbit of the Earth, the spiral increases so 
that the field is predominately azimuthal. As the result of small differences in 
the quiet-day solar wind velocity around the Sun, considerable disordering can 
be expected. 
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17. Electric field strength E . Since the electric conductivity is high, 
neither an impressed nor an induced coronal field can be maintained. The 
electric field strength for the coronal plasma may be taken as 0. 

18. The electron charge e . Its value is 
e = 1.60207 • i0 -19 coulomb . 

19. The magnetic permeability ix . The characteristic value for the 

6 

coronal plasma may be taken as die value in free space: 

\x = 4ir • 10 7 kg m coul -2 

20. Dielectric constant e. The characteristic value for the coronal 
plasma may be taken as the one in free space: 

e = — r~ • 10" 9 kg -1 m~ 3 sec 2 coul 2 
3b TT 

21. Electrical conductivity <5 . It is obtained from the equation 

3 / 

5 = 7 • 10 2 T n ohm" 1 m -1 

The values are contained in Table II. 

22. Friction coefficient a , which characterizes the effect of the 
viscous forces on the electromagnetic force. It is given by the relation 


Substituting the previously given values, the results contained in Table 
II are obtained. 

These 22 quantities may be used to characterize the coronal plasma. 
The values in Table II, which belong to T = 1. 02 • 10 s °K, may be selected 
as characteristic values. 
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CHARACTERISTIC VELOCITIES, FREQUENCIES 
AND LENGTHS 


Some velocities, frequencies, and lengths may be derived from the 22 
characteristic quantities, which are useful in obtaining nondimens ional parameters. 

1. Velocity of light: c = — - = 2.99793’ 10 8 m sec -1 

nTc ~ 

e 

2. Velocity of sound: a = \l yk T 


Substituting 

5 

7 = 3 

and the previously given values, the results contained in Table HI are obtained. 

3. Velocity of the Alfven wave: V, = 

h 

Substituting the previously given values , 

III are obtained. 

4. The plasma frequency: 

rr 

co = e / 

p n m e 

It is the frequency of oscillation due to the electric field alone. Sub- 
stituting the previously given values, the results contained in Table III are 
obtained. 


B 

the results contained in Table 


5. 


The cyclotron frequency: 



TABLE in. CHARACTERISTIC VELOCITIES, FREQUENCIES, AND LENGTHS 


a 

s 

47. 9834 

75.3509 

76. 1099 

322. 5270 

345.66(30 

684.2160 

760. 1440 

769.4130 

1290.8700 

1396.6500 

2164. 9000 

2525. 2700 

3027. 2600 

asm 

373453.00 
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It is equal to the angular frequency with which particles of mass m and 
charge e gyrate in a cyclotron. Substituting the previously given values, the 


results for the electron 
obtained. 


u) and for the proton c o contained in Table III are 


ce 


cp 


6. The electrical decay frequency: 

_J _ 6 
W e D “ e 

It prescribes the rate of electrical energy converted into joule heat. 
Substituting the previously given values, the results contained in Table III are 
obtained. 


7. The friction frequency: 


a _ a 
f mn ~ p 

It is the frequency of the oscillation if the force of friction is the only 
external force. Substituting the previously given values, the results contained 
in Table III are obtained. 


8. The Debye length: 


I kT 
D v ne 


It is a measure of the distance over which an excess electrical charge 
may be appreciably different from zero. Substituting the previously given values, 
the results contained in Table HI are obtained. 


9. The mean free path: 
___ M 

i . = 1.255 \Ty — 
f a p 


Substituting y 


5 

3 * 

M 

1.6202 — . 

ap 
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It is a measure of the distance traveled between collisions of neutral 
particles . Since the typical length of the flow field is of the order of the mean 
free path, or smaller, the collision theory should be applied. Substituting the 
previously given values, the results contained in Table in are obtained. 

10. The Larmor radius: 



It is a measure of the radius of the helical path in a charged particle 

in a magnetic field. Substituting the previously given values, the results for 

the electron l . and for the proton l „ are contained in Table III. 
f e F ip 


NONDIMENS I ONAL PARAMETERS OF CORONAL DYNAMICS 


The nondimens ional parameters which characterize the flow of a fully 
ionized plasma are obtained by dimensional analysis from the 22 physical 
quantities given earlier in this report. In the system of units adopted, there are 
five independent fundamental units: mass, time, length, electrical charge and 
temperature. Then, by the ir -theorem of dimensional analysis, 17 nondimen - 
s ional parameters may be formed. The following five fundamental units are 
suggested: 

£ = 10 9 m 

o 

p = 1.80* 10 16 kg m" 3 
o 

t s 13 330 sec 
o 

e = 1.60207 * 10” 19 coulomb 

T = 10 6o K 
o 

Nondimens ional parameters of all other quantities listed under Typical 
Physical Quantities are then derived in terms of these fundamental units. They 
will be denoted as n. . 

i 
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i. 


The nondimens ional parameter for the velocity: 



2. The nondimens ional parameter for the pressure: 



2 

3. The ratio of specific heats: 


n 3 = 7 = 


_ _J2_ 


4. The Mach number : 

u 

n, = M = - 

5. The Reynolds number: 

„ p ui 

n 5 = R = - — 

& e p 
v 

6 . The Prandtl number: 


n 6 = P r = 


p c 

v p 


7. The Froude number: 

u 2 

n 7 = f = —r 
1 r gi 

8. The Schmidt number: 


n 8 = 


A. 

P D, 
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9. The relativity parameter: 


u‘ 

n 9 = R c = V 


10 . The electric decay parameter: 


n 10 = t co = t — 

1U e e 


11. The plasma frequency parameter: 


Iln = to > = te / — 

“ p n me 


12, The frictional frequency parameter: 


a 


II 12 = t co. # t — - 
1 f mn 


13. The electrical frequency parameter: 

E 

1^13 = = - rr 

E u u H 
e 

14. The magnetic pressure number: 


n “ = R h ' "7 


ju H 
e 


15. The magnetic Reynolds number: 


II 15 = R - = u i 

o e 


ui 

v 


V = 

h 
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16. The radiation pressure parameter: 

P, 


n i6 


r 

P 


a T 4 
r r 

3p 


17. The radiation flux parameter: 


D a 

n 17 = Ti 


The following seven parameters have considerable influence on the flow 
problems of the corona: 

H 2 , y , M , R g , P r , R h , and R 6 . 

The remaining seven are unimportant parameters which may be neglected 
F , S , R 10 , 1% , n 12 » n 16 , and IIj 7 . 

r c 

Substituting the previously given values, the results for the important 
nondimens ional parameters are given in Table IV. 

TABLE IV. NONDIMENSIONAL PARAMETERS FOR THE CORONAL PLASMA 


r 

R 

© 

n 2 

n 4 = M 

*-*• 

n 6 = P r 
10“ 2 

n U - ®h 

n , 5 « R 6 
10‘ 9 

i . 58 

445.9450000 

0. 0516243 , 

3. 496270 

6.77280 

1162.580000 

1.53521 

1.99 

121.6590000 

0.0988077 

2. 358820 

6 . 77279 

380. 133000 

2 . 82200 

2.00 

146.5150000 

0.1010030 

2.329310 

6.77279 

369. 106000 

2.88472 

4.50 

5 . 2654100 

0.4747850 

1. 120710 

6 . 77282 

18.511100 

6.65736 

4.67 

4. 8592600 

0. 4942400 

1.070660 

6.77279 

16.524700 

6.79625 

6.80 

2.1643000 

0.7411640 

0.711583 

6 . 77281 

5. 570630 

7.74194 

* 7.20 

1.9487400 

0 . 7835430 

0.661667 

6.77279 

4.761170 

7. 89263 

7.25 

1.9294700 

0. 7848900 

0.651804 

6.77281 

4.741860 

7.85198 

9.56 

1.2487000 

0. 9761290 

0.445223 

6 . 77280 

2.287770 

8.46466 

10.00 

0.9945580 

1.0918600 

0 . 435902 

6.77280 

1.866610 

9 . 05209 

12.60 

0.5608700 

1.4553100 

0.310081 

6.77280 

0.901651 

10. 96430 

13. 70 

0.7605530 

1.2498500 

0.283891 

6.77280 

0. 889260 

8.92412 

15.00 

0.3261180 

1.9097200 

0.234148 

6.77280 

0.491019 

13.63560 

214. 00 

0.0578067 

4. 5341600 

0 . 989044 

6.77280 

2384. 120000 

4. 36873 
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